Dl and D2 lines in ^^Li including QED effects 
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Accurate theoretical predictions including leading QED corrections for 2^Pi/2 ~ 2^51/2 (Dl) and 2'^P-j^2 ~ 
2^<S'i/2 (D2) transition energies have been obtained for ®'^Li isotopes. Our results for the Bethe logarithms 
lnfco(2^S') — 5.178 169(4), and for mass polarization corrections Alnfco(2^S') = 0.11381(3) are in dis- 
agreement with ones obtained recently in the literature. In contrast, results lnfco(2^P) = 5.179 81(7) and 
Alnfco(2^P) = 0.1113(5) are in good agreement with them. From our theoretical predictions and recent 
measurements of '''^Li D-lines at NIST, we determine the mean square charge radius difference between ^Li 
and ®Li nuclei, in agreement with determinations based on the 3^5*1/2 — 2^51/2 transition, what demonstrates 
consistency of atomic spectroscopy determination of fundamental properties of nuclei. 



PACS numbers: 31.15.ac, 31.30.J- 
I. INTRODUCTION 

Accurate spectroscopy of the lithium atom allows one to 
test computational approaches to many electron systems in- 
cluding quantum electrodynamic (QED) effects. Compari- 
son with measured values not only verifies theoretical pre- 
dictions, but also experimental values as many discrepancies 
have been reported in the literature. A very good example is 
the recent measurement of D-lines performed at NIST ^ IH, 
which was stimulated by a long standing discrepancy among 
theoretical and experimental results for the fine structure split- 
ting 2P3/2 — 2P1/2, and its difference between isotopes. This 
has been calculated in many works and only recently, it was 
discovered a significant effect of the hyperfine mixing of P- 
levels I3J. From the experimental side, the presence of the 
nearby level, whose natural width overlaps with measured 
transitions, may contribute significantly to the lineshape and 
to the linewidth When it is taken into account in deter- 

mination of line positions, the most recent experimental result 
|[ll|2l becomes in a very good agreement with theoretical pre- 
dictions. Moreover, assuming that theoretical predictions and 
recent measurements [2J are as accurate as claimed, the mean 
square nuclear charge radii difference between ^Li and '^Li 
can be determined. Our result and that presented in experi- 
mental work [2], are in an agreement with the one obtained 
from 35* — 2S transition. This strongly supports the atomic 
spectroscopy determination of nuclear charge radii, not only 
for stable nuclei, but also for very short lived ^^Li and ^^Be 
isotopes, which can not be studied by other means. 

In this work we present the most accurate calculation of 
relativistic and leading quantum electrodynamics corrections 
to lithium D-lines. The calculational approach is similar to 
that developed previously by Drake and Yan in Ref. f?], with 
some differences regarding evaluation of three electron inte- 
grals. We use recurrence relations, derived in (Q. They are 
sufficiently stable and fast to obtain nonrelativistic energy lev- 



els to 10~^^ precision in a one day calculations ff]. Apart 
from nonrelativistic energy, other contributions to energy, rel- 
ativistic and QED ones require evaluation of more compli- 
cated matrix elements, what is probably, the most challenging 
part of these calculations. Here, we use an approach devel- 
oped in |@,@], and obtain results in good agreement with that 
by Drake and Yan in f^, with exception of Bethe logarithms 
and it's mass polarization corrections. They have not so far, 
been verified with competitive accuracy, and indeed our result 
for 25* state is not in an agreement with that obtained in [5^. 



II. EXPANSION OF THE ENERGY LEVELS 

We follow the approach used in the former calculations on 
the relativistic and QED effects for three-electron systems 
8]. Here we concentrate on Dl and D2 transition energies and 
the related isotope mass shifts. 

Energy levels are expanded in a power series of the fine 
structure constant a and the reduced electron mass to nuclear 
mass ratio ?; = —fi/niN 

E = [£(2,0) +^^(2,1) +^2^(2,2)] 

+ ma^ [£(4,0) + j-(4,o) ^ ^ (^(4,i) ^ _^(44))] 
+ ma5 [£(5,o)+^£(5,i)] (1) 
+ maS£:(6,o)^^^7^(7,o)^^^^^ ^2) 

with dimensionless spin independent (centroid) £'(™'") and 
spin dependent J7(m,n) coefficients. All neglected terms de- 
noted by ". . ." involve higher powers of a and the mass ratio 
?; and thus are small in comparison to the accuracy of the final 
results. The leading contribution £(-^'") = £9 is a solution of 
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the Shrodinger equation with the clamped nucleus 

Ho* = fo^-, 
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The trial wave function 4* is a linear combination of ij] ele- 
ments, the antisymmetrized product A of the spatial function 
(f) and the spin function xttlSll 



i>s = A[(t)s x] , 

if Pa = A[(t)PaX\ 
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X = (a(l)/3(2)-/3(l)a(2))a(3). 
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with fii being non-negative integers, Wi E R+, and the sub- 
script a = 1,2,3. The matrix element of the Hq in Eq. (|4|i 
(or of any spin independent operator) can be expressed after 
eliminating spin part x as 

imol^') = (2 0(1,2,3) + 2<^(2,1,3)-^(3,1,2) 
-</)(2, 3, 1) - ^(1, 3, 2) - <^(3, 2, l)\Ho |0'(1, 2, 3)) .(10) 

In comparison to our former calculations in Ref. [8], we 
have improved in numerical results for nonrelativistic en- 
ergy of the ground state, 

£o(25) = -7.478 060 323 910 2(2) (11) 

and here for the first excited state 

Sa{2P) = -7.410 156 532 651 4(9). (12) 

They are however, about 2 orders of magnitude less accurate, 
compared to the most recent large scale calculations by Wang 
et al. 1 10]. We note, that overall theoretical predictions are 
not limited, at present, by nonrelativistic energy, but by the 
unknown m contribution. 

In the infinite nuclear mass limit, we neglect in the non- 
relativistic Hamiltonian the mass polarization term, which is 
included later as the perturbation to T-Lq 
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The first order perturbation to the wave function ^ 
V S'i'mp is defined by 
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as a correction 



originating from the mass polarization perturbation of the 
wave function. Then, the first and the second order finite mass 
nonrelativistic corrections are obtained from the mass scaling 
and the mass polarization corrections 



£•(2,1) ^ £-(2,0) 
£(2,2) ^ 
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From the Breit-Pauli Hamiltonian, we extract the spin inde- 
pendent part 



^(4,0) ^ J2 
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leading to the coefficient E^'^^'^^ — {Ti.^^'^^). The finite nu- 
clear mass correction consists of the mass scaling term 
£ms^\ the mass polarization perturbation ('H''''°')mp, and the 
explicit recoil term Tiff'''^^ coming from the Breit interaction 
between electrons and the nucleus, 

5(4.1) ^ 44,1)^ ^^(4,0)^^^ ^^(4,1)^^ (18) 

- -4E|(p')+3[^E(^'(-'^)) (19) 
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The spin-orbit coupling is given by 



(4,0) _ 
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(21) 



E ^7~3" ^'^ [9 rab X Pb - (g - 1) fab X Pa] 
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(22) 
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_^(4,l) ^ 3 ^^(4,0)^ ^ ^^(4,0)^^^ ^ ^^(4 0)^ _ (24) 



with spin matrices represented using Pauli matrices Sa = 
(Ta/2. The matrix element of spin-orbit terms are obtained 
for nPi/2 state by using 



(V^alE^^-'^-l^'') = i(<?a(l,2,3) 

c=l 

-2Q3 X [^b(l,2,3) + 0fc(2,l,3)] 

+ (Qi -Q2 + Qs) X [M2, 3, 1) + 06(3, 2, 1)] 

+ {-Qi + Q2 + Q3) X [0b(l, 3, 2) + 0fc(3, 1, 2r)] 

and forthe nP3/2 state J"(P3/2) = -1/2 J'(Pi/2). 



(25) 
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Leading QED correction of order m is given by 
'19 
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where the Bethe logarithm In /cq has the form 

In fco = 

(Eg Pa (Ho - Eq) In [2 (Ho - Ep)] Y.bPb) 

and the Araki-Sucher term 
1 



(26) 



(27) 
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lim / d'^rd>*(r) 
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Q{r - a) +ATTS^{r) (7 + Ina) 



(28) 



The finite nucleus mass correction to the leading QED con- 
tribution has the form 



c(54) _ P(5,l) i r(5a) I c(5,l) 



(29) 



The mass scaUng part is given by 

4i 
"3 
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Apart from the scaling factor three for all operators in Eq. 
(|26] |. there are additional two terms from the mass scaling of 
the Bethe logarithm and P(l/r'^) operator The direct recoil 



,^(5,1) 
C(5.1) 



term Si; ' is known as the Salpeter correction 
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and £mp^^ is due to the mass polarization correction to expec- 
tation value of the operators in Eq. (|26] |. 

The relativistic m corrections for few-electron atoms are 
very difficult to calculate, and closed formula has been ob- 
tained only for two-electron systems. The result for three- 
and more-electron systems lUlll contains divergences, which 
elimination has not yet been performed. For three-electron 
systems we use an approximate formula on the basis of hydro- 
genic values including dominating electron-nucleus one-loop 
radiative correction f 12l| 
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We neglect electron-electron radiative corrections and the 
purely relativistic corrections, as we expect them to be rela- 
tively small, of order 10%. This is the leading source of uncer- 
tainty in the theoretical predictions for transition frequencies. 
Similarly, the finite nuclear mass correction, which is signifi- 
cant for the isotope mass shift, is also estimated on the basis 
of the known hydrogenic values. Apart from the mass scaling 
and the mass polarization corrections to the Eq. (l32t , the ra- 
diative recoil and pure recoil corrections have been included 
as folows 
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Relativistic and QED contributions to the fine structure have 
been described in details in Ref. [3], and we note that the 
hyperfine interaction leads to a significant shift in the ^Li fine 
structure. 

Due to the numerical significance for transition energies, 
one estimates the m cP contribution on the basis of formulas 
for hydrogenic systems II12I1 . 



(7,0) 



(n) 



m-^^ [A6o(n)+ln(Za)-2Ai(n) 

TT TV' 

+ ln^ 2^g2 +m - ^ 
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-ml- ^ — ^ C40. 



(34) 



It involves one-, two-, and three-loop corrections, and values 
of A, B, and C coefficients may be found in {vi\. Following 
Ref. fist] these hydrogenic values of order m are extrapo- 
lated to lithium, according to 

^ [2£C^){\S,Z) + £('^\nX,Z-2)\ 
(,53(ri)+,53(r2)+53(r3))Li 



2{5^{r)) 



1S,Z 



(35) 
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for X = S, and for states with higher angular momenta 
£^'^\nX,Z) is neglected. We assume this approximate for- 
mula to be accurate to 25%. 

Beyond QED, there are corrections due to the finite nuclear 
size. The leading order correction ma'^ is given by 



(4,0) 



(36) 



where is the averaged square of the charge radius and X 
is the electron Compton wavelength divided by 2 tt. We ad- 
ditionally include a logarithmic relativistic correction to the 
wave function at the origin 



f(6.0) 
'fs,loK 



-(Z 01) ln(Z a TO Tc) 



(4,0) 



(37) 



and the finite mass correction £, 



(4,1) 



which consists of the 



mass scaling and mass polarization correction to the Eq. 
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One determines nuclear charge radii from the difference be- 
tween experimental and theoretical isotope shift by using 



exp 

Cab 



CABir^A-r'cB), (38) 
J' • (39) 



'cA 'cB 

the described finite nuclar size f ts for isotopes A and B. 



3/2Z fl6"]. It helps in jud ging on the completeness of the 
intermediates states and on the estimation of uncertainties. 

Calculations for the 2P state are more sophisticated as com- 
pared to the ground state. The dipole operator P couples the 
wave function of the 2P state to three classes of intermediate 
states S, P'^, and D which differ by the orbital angular mo- 
mentum. This fact allows to split the optimization of (5^in 
into three parts, involving 05, and 



III. BETHE LOGARITHM AND ITS MASS 
POLARIZATION CORRECTION 

Bethe logarithm is the most demanding term in accurate nu- 
merical evaluation among all operators in Eq. (|2). In this work 
we use the integral representation introduced by Schwartz U41 
with compact set of formulas jTsI] given by 



In fco = 


DJo 


/O - f2t^ 


(40) 


/w = 






(41) 
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\/l + 2w 
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fo = 




-2D, 


(44) 



The integrand in Eq. ( l40l l is a regular function with a reason- 
able good numerical convergence calculated at 100 equally 
spaced points in t e [0, 1] variable. The critical region are 
points at low t. Here, following Schwartz [14], we perform 
the expansion in small t 

fit) = fo + f2t^ + fzt^ + fAtHn{t) + 0{t^), 
h = 8ZD, U^ieZ^D. (45) 

The value of the integrand in Eq. (l40l i at t = is equal to 
/a and this helps to judge the numerical precision achieved at 
small t-region. In order to perform integration, we use poly- 
nomial interpolation of the integrand for t > 0.1. For the 
remaining region, t e [0, 0.1] we extend the expansion of the 
remainder denoted by 0{t^) in Eq. ( |45] l to higher powers of 
t, next we fit the expansion terms and finally perform the in- 
tegration. 

The function f{t) can be represented as the matrix element 
with vj/ and a pseudo-state function S'i'in given by 

fit) = -Lu{^\P\S^i,,) (46) 
(£o - Ho - 5*in = P^. (47) 

In the calculations of the Bethe logarithm for S'-states, 6'i'iu is 
of the form of a P-function in Eq. (|6]l. This function is calcu- 
lated by optimizing f{t) against nonlinear parameters of (5^i„ 
for each t individually. The special case is f{t) /uj in the limit 
a; = 0, which is established by the Thomas-Reiche-Kuhn 
sum rule for dipole oscillator strengths {P (Ho — £q)^^ P) — 



6S 



<^tjkrlii (t>s , 



■ra-n] (t)s 



(48) 
(49) 



respectively. In order to be able to use the optimization against 
nonlinear parameters, the negative contribution coming from 
25 to (J^tin is represented as a fixed sector with parameters 
taken from the ground state wave function. All the other non- 
linear parameters are obtained by optimization of the corre- 
sponding parts of fit). 

The mass coefficient S^p ^ in Eq. ( |29l ). involve the mass 
polarization correction to the Bethe logarithm. The frame- 
work of such perturbative calculations have been presented on 
the recoil corrections to the helium atom [17]. Following the 
Eq. ( l40b we can build the integral representation in the form 

A(lnfco)mp = ^"^i^^lpW-/0.mp-/2,mpt' 



D 



/^p(i) = ~u{P- ^, p) 



(50) 
(51) 



mp 



P^—TT— {J2P--Pb-iP-- Pb)) (52) 



^0 — Ho — UJ 
1 



a<b 



-P 



^0 — 'Ho — u! 
D^p = 2^Z ^(J3(r„))^p, 



(53) 



fo. 



mp 



/mp. 



/2.mp = -2 Z?„,p , (54) 



Apart from the correction to the wave function Eq. (Il4t . the 
mass polarization correction Eq. (ITjt in the denominator of 
the f{t) in the Eq. (l4Tl i corresponds to the second term of 
Eq. ( fSTT i. All terms in the above representation are evaluated 
as the mean values with numerically determined functions 



S'i'mp and (J^in. 



IV. RESULTS 

Table H] presents numerical values of the Bethe logarithms 
for the lowest doublet S and P states of lithium. In order to 
control the numerical uncertainty, we performed calculations 
with several basis sets successively increasing the shell pa- 
rameter (see eg. [7]). Our result for In fco and A In ko agrees 
well with former ones based on the integral representation 
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TABLE I: Bethe logarithm and the mass polarization correction in 
the lithium atom 



Shell 


In ko 


A In fco 


25* 






7 


5.178 067 3 


0.113813 


8 


5.178 145 2 


0.113 776 


9 


5.178 1618 


0.113 835 


10 


5.178 167 9 


0.113 807 


oo 


5.178 169(4) 


0.11381(3) 


Pachucki, Komasa (2003) [15] 


5.17817(3) 


0.114(3) 


Yan et al. (2008) [5J 


5.178 28(1) 


0.113 05(5) 


2P 






7 


5.178 974 


0.11609 


8 


5.179 549 


0.11316 


9 


5.179 635 


0.11196 


10 


5.179 780 


0.11148 


oo 


5.179 81(7) 


0.1113(5) 


Yan et al. (2008) [5] 


5.179 79(6) 


0.1112(5) 



with ECG functions fT3\, but are more accurate. An incon- 
sistency is observed with the resuh obtained in Hylleraas ba- 



sis by Yan et al. f3\. These authors use a discrete variational 
representation of the continuum in terms of pseudo states to 
cover a huge range of distance scales [fisl]. At the present 
level of theoretical predictions, the resulting differences of 15 
MHz for transition energies 2P — 2S and 1 kHz for isotope 
mass shifts are much smaller compared to the uncertainties 
from estimations of higher order QED corrections. For 2P 
state our results confirm those by Yan et al. [5^ and both have 
comparable uncertainties. 

Table HI] presents results for various spin independent 
g{™.,n) (jimensionless coefficients. Except for ma^, ma^ and 
ma^rj terms, all uncertainties have been obtained from the 
analysis of the numerical convergence as a function of the dif- 
ferent basis size for the wave function, since operator struc- 
ture at given order in a and 77 are included in a complete way. 
Complete result for £(^'0) is not yet known, it is estimated by 
what we think, the dominating one-loop radiative correction, 
see Eq. (|32] l. The uncertainty of £(6 '^) is assumed to be 10%, 
what is a dominant source of the overall uncertainty. Anal- 
ogous estimations of and f^^^^^ are established at the 
level of 25%. Next columns of Tablellllpresent (centroid) tran- 
sition energies vq, 1/7, and the (centroid) isotope shift Ai^g?- 
One observes that each power of a or 77 gives significantly 
smaller contributions, so the expansion in a and -q is physi- 
cally well meaningful. 



TABLE IL Spin independent contributions (centroid) to the transition energy v{2P — 2S) for ®Li, '^Li and to the isotope shift. 



energy 


^(m,n) 


!^6(MHz) 


!/7(MHz) 


Al'67(MHz) 


mop 


0.067 903 791259 0(16) 


446 785 483.5(1) 


446 785 483.5(1) 




4 

ma 


0.267 6121(4) 


93 765.1(2) 


93 765.1(2) 




ma^ 


-3.469(3) 


-8 850.(8) 


-8 870.(8) 




mo^ 


-14.4(1.4) 


-269.(26) 


-269.(26) 




ma' 


217.(54.) 


30.(7) 


30.(7) 




ma^rj 


0.123 007 926 0(3) 


-73 826.6 


-63 293.1 


-10 533.510 5 


ma^r) 


-0.267 591(2) 


8.6 


7.3 


1.220 2 


ma^r) 


1.134(3) 


-0.3 


-0.2 


-0.037 7(1) 


ma^ri 


-46.(12) 






0.011(3) 


2 2 
ma T] 


-0.004 870(4) 






-0.070 7 


f(i.O) 


-1.045608 r? 


-16.0 


-14.0 




C67 








-2.465 8 


Total 




446 796 306.(28) 


446 806 840.(28) 


-10 532.388(3) 



Table Hnlpresents results for contributions to the fine struc- 
ture splitting and the corresponding isotope shift. The result 
for mixing of 2P3/2 and 2Pi/2 due to hyperfine interaction 
comes from our former paper [3|. Other contributions can in 
principle be deduced from the individual operators given in 
Ref. 0,13], but results presented here are recalculated with 
higher numerical precision, and they are the most accurate 
among ones available in the literature (except ma^ [10]). We 
note that the present numerical precision for expansion coef- 



ficients is high enough, so that the uncertainty comes exclu- 
sively from neglected higher order corrections. These correc- 
tions are not expected to be very small, therefore in compari- 
son to experimental results by Sansonetti et al. theoret- 
ical predictions for the fine structure are much less accurate. 
However, the isotope shift in the fine splitting is not so much 
sensitive to higher order terms, as they cancel out in the dif- 
ference. As it was pointed out by Yan et al [5j, various ex- 
periments and theoretical predictions are all in disagreement 
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between themselves. Only recently, after correcting for the start to agree with very recent experimental value obtained by 
hyperfine mixing the theoretical predictions of —0.544 7(1) NIST group -0.531(24) ll. 



TABLE IIL The fine structure in ®Li, ^Li and tlie fine structure isotope shift '"'^Li in MHz. The Total line includes only numerical uncertainties, 
and not uncertainties from the unknown m a'' correction. This correction is negligible only for the isotope shift in the fine structure. 



energy 



-3/2.F('"'")(2Pi/2) t.6(2P3/2 - 2Pi/2)(MHz) ;/7(2P3/2 - 2Pi/2)(MHz) Ai.67(2P3/2 " 2Pi/2)(MHz) 



ma 

mixing [J] 
Total 



0.028 693 979(3) 
0.08718(2) 



10 053.712 6(11) 
-2.7861(5) 
0.012 17 
10 050.938 6(12) 



10053.712 6(11) 
-2.388 6(4) 
0.159 16 
10 051.483 2(12) 



-0.397 5(6) 
-0.146 99 
-0.544 7(1) 



Experiment |T] 
Experiment [2] 



10 052.799(22) 
10 052.779(17) 



10 053.393(21) 
10 053.310(17) 



-0.594(30) 
-0.531(24) 



Transition energies for Dl and D2 lines in Table llVlare ob- 
tained from the data presented in Tables HH and Hill As for the 
fine structure, theoretical predictions are much less accurate 
in comparison to experimental values 01 lid due to unknown 
complete m contribution. Nevertheless, our results are in 



good agreement with previous theoretical calculations by Yan 
and Drake in The small difference comes mainly from 
previously mentioned difference in ground state Bethe loga- 
rithms in Table |T] 



TABLE IV: Transition energies for Dl and D2 lines. Comparison of theoretical and experimental values. 



energy 



2Pi/2 -25'i/2(MHz) 



2P3/2 -25'i/2(MHz 



This work 
Experiment [5] 
Experiment |T] 
Experiment \2} 

This work 
Theory [5] 
Experiment |_5] 
Experiment |T] 
Experiment [2] 



446 789 589.(20) 
446 789 596.091(7) 



446 800123.(20) 
446 800 142.(30) 
446 800 130.61(42) 
446 800129.853(6) 



«Li 



^Li 



446 799 640.(20) 

446 799 648.927(21) 
446 799 648.870(15) 

446 810175.(20) 
446 810 193.(30) 
446 810 189.18(42) 
446 810183.289(20) 
446 810183.163(16) 



Table IVl provides a new determination of the difference in 
squared charge radii Sr^ from a new measurement of the iso- 
tope shift of the Dl and D2 hnes for ^Li and *^Li tM- The final 
results 0.705(3) from Dl (the most precise in the literature) 
and 0.700(9) from D2 lines are in a very good agreement. 



They both agree with the determination based on 3S' — 25* 
transition 0.731(22) |20], but not so perfectly. The difference 
is about one a. It nevertheless stands as a confirmation of 
atomic spectroscopy approach to determination of the nuclear 
charge radii. 



V. SUMMARY 

We have calculated lithium D-lines including leading QED 
effects with high numerical precision. All results, but Bethe 



logarithms are in agreement with previous ones obtained in 
the literature. A discrepancy is observed for the Bethe log- 
arithm of ground state, which needs further verification. Al- 
though all these Bethe logarithms are somehow close to hy- 
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TABLE V: Summary of the isotope shift determination of the "Li charge radii from Dl and D2 lines with respect to ^Li, 
r(''Li) = 2.39(3) fm UJ], the first uncertainty of Uthc comes from unknown higher order terms, the second uncertainty is due 
to the atomic mass. Ce? = —2.465 8 [MHz fm"^] 





D-line 


!/exp[MHz] Ref. [2] 


t'tho[MHz] 




rc[fm] 


2Pl/2 


— 251/2 


-10 533.763(9) 


-10 532.023 7(28)(2) 


0.705(3) 


2.533(28) 


2^3/2 


— 2Si/2 


-10 534.293(22) 


-10 532.568 2(28)(2) 


0.700(9) 


2.533(28) 


3S - 


25* [20] 






0.731(22) 


2.538(28) 



drogenic value, there is no any computational approach which 
takes advantage of it. One always has to represent the interme- 
diate states very accurately and this is a numerically demand- 
ing task. All the other corrections can be obtained with much 
higher accuracy at almost no computational costs. In compar- 
ison to experimental D-lines, agreement is observed but theo- 
retical predictions are anyway much less accurate. The main 
reason, which is a principal problem in improving theoretical 
results for lithium, is the complete calculation of m cor- 
rections, which so far has been performed only for He ||2T]| . 
This is a challenging, but what we think, a feasible task, with 
computational approach based on explicitly correlated basis 
functions, such as Hylleraas, Slaters or Gaussians. This cor- 
rection cancels out in the isotope shift in fine structure, and 
here theoretical results are in agreement with the most recent 
result of Sansonetti et al |[ll|3l. An additional results of this 
work is the determination of ^'^Li nuclear charge radius dif- 



ference from the measured and calculated isotope shift. Our 
result is maybe not in perfect, but in one cr agreement with that 
obtained from 35 — 25 transitions, what clearly demonstrates 
the applicability of atomic spectroscopy methods for accurate 
determinations of nuclear charge radii. We should mention, 
nevertheless, that for the simplest possible atoms, the hydro- 
gen and the muonic hydrogen, a significant and unresolved 
discrepancy is observed for the proton charge radius [22J . 
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